Universal bosonic tetramers of dimer-atom-atom structure 
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Unstable four-boson states having an approximate dimer-atom-atom structure are studied using 
momentum-space integral equations for the four-particle transition operators. For a given Efimov 
trimer the universal properties of the lowest associated tetramer are determined. The impact of 
this tetramer on the atom-trimer and dimer-dimer collisions is analyzed. The reliability of the 
three-body dimer-atom-atom model is studied. 

PACS numbers: 34.50.-s, 31.15.ac 



INTRODUCTION 



It took more than 30 years till Efimov's prediction 
for the existence of weakly bound three-particle states 
with asymptotic discrete scaling symmetry [l| was con- 
firmed in the cold-atom physics experiments [2|- This 
discovery raised the interest in few-body systems with 
resonant short-range interactions, both experimentally 
and theoretically. After the three-body system where 
semi- analytical results have been obtained as summa- 
rized in Ref. [3|, the next step in the complexity to ex- 
plore the Efimov physics is the system of four identical 
bosons where several numerical studies |4l4l0| are already 
available. Although there is no four-boson Efimov ef- 
fect [ll|, HJI , the properties of the four-boson system are 
strongly affected by the three-boson Efimov effect. The 
observables of the four-boson collisions show the same 
discrete scaling symmetry and are related to the respec- 
tive trimer binding energies in a universal way provided 
the size of the involved few-body bound states greatly 
exceeds the interaction range [a, Q- Furthermore, be- 
low each Efimov trimer in a certain regime there are two 
tetramers [5|, |6| . The tetramers associated with the ex- 
cited trimers are unstable bound states (UBS) (iG. Ha. [l4| 
and lead to resonant effects in the four-boson reactions. 
First numerical calculations of the four-boson recombi- 
nation ^j and dimer-dimer collisions |15i] in the adia- 
batic hyperspherical representation were followed by the 
work of Refs. [3, [la, [I^l employing the integral equa- 
tions for the transition operators that were solved in the 
momentum-space framework. In the latter case the uni- 
versal limits for the scattering observables and tetramer 
properties such as widths and intersections with thresh- 
olds were obtained with considerably higher accuracy; 
furthermore, remarkable resonant effects were predicted 
in the atom-trimer collisions as well [IJ] . The theoretical 
results [a, fl5l - [l7{ are roughly consistent with the existing 
data from the dimer-dimer and four-atom recombination 
experiments |18l - [2l| performed in a regime that is not 
strictly universal. 

In the present work we will study a different type of 
bosonic tetramers, namely, the ones of the approximate 
dimer-atom-atom structure. Their existence was pre- 
dicted in Refs. 0, [l^ as a consequence of the three-body 



Efimov effect in the three-body system made off a dimer 
and two atoms. This is illustrated in Fig. [T] where we 
schematically show the four-boson energy {E) spectrum 
as function of 1/a with a being the two-boson scattering 
length. We consider the regime around the special value 
a = a"^ that corresponds to the nth Efimov trimer being 
at the atom-dimer threshold, i.e., hd = &« where hd {b„) is 
the dimer (nth trimer) binding energy relative to the free 
particle threshold. If the trimer is a true bound state with 
zero width, i.e., if deeper dimers are absent (which is not 
the case in typical experiments |18l42l| ). the atom-dimer 
scattering length Ad is infinite at a = a„. Sufficiently 
close to _E = — 6„ and a = af^ where Ad greatly exceeds 
the dimer size being of the order of a, one may expect 
to mimic some properties of the few-boson systems using 
a model that considers the dimer as a pointlike particle. 
Then in the effective three-body system consisting of a 
dimer and two atoms there are two atom-dimer pairs with 
infinite two-body scattering length Ad- In such a three- 
body model of the four-boson system the three-body Efi- 
mov effect occurs, however, with a very large discrete 
scaling factor e'^/^o w 2.016 x 10^ [l]. The resuhing Efi- 
mov states accumulate at E — —hd- Their dimer-atom- 
atom structure is only an approximation but no attempt 
has been made so far to describe them rigorously as four- 
boson states. This will be done in the present work us- 
ing exact four-particle equations. We label the tetramers 
of this type by two integers (n, m) where n refers to 
the associated Efimov trimer and m = 0, 1, 2, . . . distin- 
guishes between different states existing around a — af^ 
with fixed n. The (n, ni)th tetramer intersects the nth 
atom-trimer threshold at a = a*^ ^ and the dimer-atom- 
atom threshold at a = aj^^. The calculation of these 
tetramers is technically very difficult task owing their 
weak binding and very large e'^'^"; our results will be 
limited to m = states. Furthermore, all these tetramers 
lie above the dimer-dimer threshold and above all lower 
atom-trimer thresholds n' with n' < n. Therefore the 
considered tetramers are UBS with finite width. We will 
extract their properties using rigorous four-particle scat- 
tering calculations. Performing in addition the three- 
body calculations we will establish the limitations of the 
dimer-atom-atom model. 

In Sec. m] we shortly recall the technical framework. 
In Sec, mil we present results for tetramer properties and 
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FIG. 1. (Color online) Schematic representation of the four- 
boson energy spectrum in the vicinity of the dimer-atom-atom 
(dotted line) and the nth atom-trimer (dashed line) threshold 
intersection at a = af, . The tetramers (solid lines) accumu- 
late at a = a^ and E — —bd but for simplicity only two 
states (m — 0, 1) are displayed; only qualitative relations are 
preserved. The (n, rn,)th tetramer intersects the dimer-atom- 
atom and the nth atom-trimer thresholds at a = an,m ^nd 
tf = in,m) respectively. 



their effect on the atom-trimer and dimer-dimer scatter- 
ing observables. We summarize in Sec. IIVI 



II. FOUR-BOSON SCATTERING EQUATIONS 

To describe the colhsions in the fom'-boson system we 
use exact Alt, Grassberger and Sandhas (AGS) equations 
[22| for the transition operators U^a in the symmetrized 
form 

hill — P34{GotGo) + P34UiGotGoUii + U2GotGol^21, 



U21 = (1 + P34)(GoiGo)-l + (1 + P34)UiGotGoUii, 

Ui2 = (GotGo)"' + P34UiGotGoUi2 + C/aGoiGoWa: 



U22 = (l + P34)UiGotGoUi2. 



(la) 

(lb) 

'(Ic) 
(Id) 



Here Go is the free resolvent, P34 is the permutation op- 
erator of particles 3 and 4, t = v+vGot is the two-particle 
transition matrix derived from the potential v, and Ua 
are the symmetrized AGS operators [23] for 3 -f 1 (a = 1) 
and 2-1-2 (a = 2) subsystems. We solve the AGS equa- 
tions ll]) in the momentum-space framework using the 
partial-wave decomposition. In the present study only 
the states with zero total angular momentum and posi- 
tive parity (0^) corresponding to the quantum numbers 
of the tetramers need to be considered; among them the 
states with the nonconserved angular momentum of the 
three-boson subsystem J < 3 have to be included for the 
convergence. The atom-atom interaction v is taken over 
from Ref . Q ; it is given by a separable S'-wave potential 
with one- or two-term gaussian form factors. Further 
technical details can be found in Refs. [l7l. |23|. 12411 . 



The scattering amplitudes for all elastic and inelastic 
two-cluster reactions arc obtained as the on-shell matrix 
elements of the AGS transition operators [23]; the am- 
plitudes for breakup and recombination [l7[ are given by 
the integrals involving Upa ■ 

The unstable tetramers manifest themselves as poles 
of the operators W^q in the complex energy plane at 
-Btrn - i^trnf^ with -B^^ being the energy of the 
(n, TO)th tetramer relative to the four-boson breakup 
threshold and FJ^ „ being its width. The UBS pole in 
the complex energy plane is located in one of the un- 
physical sheets that is adjacent to the physical sheet [IJ] 
and affects the physical observables leading to resonant 
effects in the four-boson scattering. Thus, the proper- 
ties of the unstable tetramers can be extracted from the 
behavior oiUpa as described in Ref. [14]. 

In the three-body dimer-atom-atom (3BDAA) model 
the corresponding reactions are described by solving the 
AGS three-body scattering equations [2^. These cal- 
culations need two pair potentials. The atom-atom 
potential Vaa = w is the same as in full four-boson 
calculations. The atom-dimer potential Vad is taken 
in the momentum-space representation as (k'|tiad|k) = 
exp(— /c' /A^) Adexp(— fc^/A^). Its strength Xd and mo- 
mentum cutoff (inverse range) parameter A^ for each a 
are adjusted to reproduce Ad, (fin — bd) and atom-dimer 
effective range predicted by the original three-boson cal- 
culations using Vaa- As a consequence, Vad yields an ac- 
curate description of the low-energy atom-dimer scatter- 
ing. We note that the range ofvad with l/A^ « a greatly 
exceeds the range of Vaa that is much smaller than a. 



III. RESULTS 

Unless stated otherwise, our results refer to rigorous 
four-boson calculations. We present them as dimension- 
less ratios that are independent of the short-range inter- 
action details in the universal limit. For this one needs 
to consider reactions involving highly excited Efimov 
trimers where the finite-range effects become neg ligible. 
As shown in previous calculations [3, M, [3 [la, [l3 , '^ ^ 3 
is sufficient for high accuracy (note that our nomencla- 
ture starts with n — for the ground state). This is 
fully consistent with our present results. For example, for 
the (n, 0)th tetramer intersection with the corresponding 
atom-trimer threshold we obtain 030/03 w 1.609 and 
o^o/af « 1.608 using one-term form factor in the po- 
tential V and flg o/flg « 1.607 and a| Q/af « 1.608 using 
two-term form factor, respectively. Thus, we conclude 
that for large n the universal ratio is 



alJat = 1.608(1), 



(2) 



corresponding to Ad/al^Q — 5.491(3). The calculation 
of the tetramer intersection with the dimer-atom-atom 



threshold, 



-'n,0' 



turns out to be numerically more com- 



plicated and less accurate due to its proximity to aji^ . We 
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FIG. 2. (Color online) Position of the (n, 0)tli tetramer 
relative to the nth atom-trimer threshold (top) and its width 
(bottom) as functions of the two-boson scattering length a 
and the atom-dimer scattering length Ad. The results of the 
full four-body (solid curves) and the 3BDAA (dashed-dotted 
curves) models are compared. The thin dotted line represents 
the zero value. 



have found that afj o/'^n > 0.9999. The higher tetramers 
(n, 77T, > 1) may exist only at extremely large atom-dimer 
scattering length Ad, i.e., extremely close to a = ajj, and 
cannot be resolved in our four-boson calculations. Thus, 
with increasing vn the ratios a^m>i/a$^ i'^i.m/'^n) very 
rapidly approach 1 from above (below). For example, our 
estimations using the 3BDAA model are 

l-alJat^l.2il)xlO-', (3) 

corresponding to Ad/af^Q = —1.8(1) x 10^, and 

an,i/«n - 1 = 1-3(1) X 10-^ (4) 



corresponding to A^/a*, 



1.7(1) X 10*^ and (6„ 



ba)/ba = 2.6(2) x IQ-i^. 

The universal values ^ and ([3]) indicate that the 
(n,0)th tetramer as UBS exists at 0.6219 < a^/a < 
1.0000012. We show in Fig. H the a-evolution of the 
(n, 0)th tetramer properties, i.e., its relative distance to 
the nth atom-trimer threshold {bn — B!^ o)/bn and its 
width rfj q/26„ as dimensionless universal quantities. We 
include also the predictions of the 3BDAA model and 
demonstrate that it is reliable only for the tetramer po- 
sition at a'^/a > 0.96 where |Ad|/a > 50. Although the 
3BDAA model describes well the nth trimer binding and 
the low-energy atom-dimer scattering, it does not sup- 
port lower atom-trimer channels and it is obviously not 
appropriate in the dimer-dimer channel since the 3BDAA 



model treats the two dimers asymmetrically. However, 
these channels are decisive for the tetramer width since 
in their absence FJ^ g vanishes. Thus, the failure of the 
3BDAA model for FJ^ g is not surprising. Furthermore, 
the pointlike dimer approximation may only be reason- 
able when its size, being of order of a, becomes negligible 
as compared to \Ad\ (that determines also the size of the 
shallowest trimer if Ad > 0), i.e., for \Ad\ >> a. Un- 
der this condition the 3BDAA model may provide a rea- 
sonable approximation of the full four-boson model for 
particular observables like the atom-trimer elastic scat- 
tering cross section and the tetramer position (6„ — Bj'J g), 
but not for reactions involving the dimer-dimer or lower 
atom-trimer channels. The intersection point estimations 
in Eqs. ^ and ^ refer to |^d|/a > 10'^ and therefore 
should be reliable whereas the 3BDAA model fails for 
^n,o/^n ^s can be seen in in Fig. [3 

The tetramer lies very close to the atom-trimer thresh- 
old, even the largest deviation around af^/a = 0.8 is only 
{bn — B^ oi/bn ~ 2.26 x 10"'*. For comparison, the dif- 
ference between the atom-trimer and dimer-atom-atom 
thresholds at this point is (5„ — bd)/bn ~ 8.86 x 10~^. 



Note that in a narrow interval around a — 



*n,0 



the 



tetramer is slightly above the atom-trimer threshold, i.e., 
B^fi < bn- The width FJ^ q vanishes at a = a^ g while 
at a > a^ Q the (n, 0)th tetramer becomes an inelastic 
virtual state (IVS) ;fl|. Negative width F^ g < in the 
TVS case implies the change of the energy sheet. The 
IVS corresponds to the pole of the transition operators 
Upa in the complex energy plane on one of the nonphys- 
ical sheets that is, unlike the one of UBS, more distant 
from the physical sheet [ij]- The IVS has a visible im- 
pact on the physical collision observables only when it is 
located very near to the scattering threshold [l3| ; an ex- 
ample referring to the four-boson system can be found 
in Ref. [lj|. In contrast, the UBS always leads to a 
resonant behavior. As example we show in Fig. |3] the 
S'-wave phase shift 63 and inelasticity parameter rjs for 
the dimer-dimer and atom-trimer scattering at the total 
four-boson energy E « —B'^q and af^/a — 0.8. It is in- 
teresting to note that only the dimer-dimer phase shift 
increases by 180 deg while the inelasticity parameter ex- 
hibits a rapid variation also for the collisions of atoms and 
trimers in the {n — l)th Efimov state. Thus, in these two 
cases the inelastic reactions are significantly enhanced by 
the tetramer UBS whereas a pronounced resonant peak 
in the elastic cross section is present in the dimer-dimer 
channel only. The (71, 0)th tetramer has very little impact 
for the collisions of atoms and trimers in the (n — 2)th 
(and lower) Efimov state. 

The (71, 0)th tetramer intersection with the 71th atom- 
trimer threshold manifests itself most prominently in the 
atom collisions with the nth trimer at vanishing relative 
kinetic energy. The corresponding atom-trimer scatter- 
ing length An is shown in Fig. [1] as a function of the 
two-boson scattering length a. An exhibits a resonant 



behavior with |Ini^„| having a peak at a 
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FIG. 3. (Color online) S-wave phase shift and inelasticity 
parameter for the dimer-dimer (solid curves) and atom-trimer 
scattering where the trimer is in the in— l)th (dashed-dotted 
curves) or in the (n — 2)th (dotted curves) Efimov state. The 
energy regime around the (n, Q)th tetramer state at a^ja = 
0.8 is shown. 
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FIG. 4. (Color online) Real (solid curve) and imaginary 
(dashed-dotted curve) parts of the atom-trimer scattering 
length An as functions of the two-boson scattering length a 
and the atom-dimer scattering length A^. Results for ReA^ 
obtained using the 3BDAA model are shown as dashed curve. 
The thin dotted line represents the zero value. 
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FIG. 5. (Color online) Inelastic reaction rates for vanishing 
relative energy atom scattering from the nth trimer as func- 
tions of the two-boson scattering length a or the atom-dimer 
scattering length Ad- Rates for the dimer production (solid 
curve) and the trimer relaxation to the (n — l)th (dashed- 
dotted curve) and to the (n — 2)th (dotted curve) Efimov 
state are compared. 



increase of ReA„ shown in Fig. linear a'^^/a — 0.99 is 
not related to the tetramers; it is due to the scaling of 
the elastic cross section with the spatial size of the trimer 
that grows with increasing Ad and decreasing (6„ — 6^), 
i.e., KeAn ^ Ad ^ {bn — bd)~^^^. This is evident in 
the inset of Fig. U] where the ratio KeAn /Ad exhibits 
no resonant behavior near af^/a = 0.99. The inelastic 
reactions quantified by ImAn are not significantly en- 
hanced by the increased size of the trimer. The failure 
of the 3BDAA model, especially for Im An and the res- 
onant peak position, can be expected based on Fig. [2] 
and its discussion. We show in Fig. U the 3BDAA pre- 
diction only for ReA„ at larger Ad/ a values. Only at 
Ad/a > 50, i.e., a^/a > 0.96, it agrees with the four- 
boson result within 7% or better. In the same regime the 
3BDAA model underpredicts Im An almost by a factor of 
40 as compared to the four-boson result. As estimated in 
Eq. (|4]) using the 3BDAA model, the An resonance due 
to the (n, l)th tetramer is expected to take place much 



closer to a 



Finally we note that rapid variations of 



An near af^/a — 0.15 are caused by the proximity of the 
dimer-dimer threshold and the Efimov tetramer (of a dif- 
ferent type) that intersect the nth atom-trimer threshold 
at ai/a = 0.14730(1) [3 and a^/a = 0.14706(1) [g, re- 
spectively. The An evolution at Un/a < 0.15 is presented 
inRef. [14]. 

The inelastic reaction rate for the atom scattering 
from the nth trimer in the ultracold limit is given [9| 
by /?„ = — (167r/i/3m)Imyl„ where m is the boson mass. 
/3„ has contributions corresponding to the dimer produc- 
tion /3'/f- and the trimer relaxation /3" to more deeply 



bound states n', i.e. 



Pn = /3f + E /^»' 



(5) 



n'=0 



The respective dimensionless quantities Pffm/ha'^ and 
/3^ m/ha'^ are compared in Fig. [5] Consistently with 
our previous findings for two-cluster collisions [3, [ig] 
and four-atom recombination [17], in most of the con- 
sidered regime the dominating inelastic reaction is the 
one leading to the most weakly bound final channel, i.e., 
the dimer-dimer channel in the present case. However, 
in a narrow interval around a^/a = 0.38 an interest- 
ing and unusual situation takes place: the dimer pro- 
duction rate /S^'' is strongly suppressed and the trimer 
relaxation process becomes the most important inelas- 
tic reaction. Another exceptional situation can be seen 
around af^/a = 0.52 where /3"~^ has a minimum such 
that /3^-i < /3;'-2. 

The considered regime is not yet explored experimen- 
tally. According to our predictions, the (n, 0)th tetramer 
most clearly could be seen as a dimer production rate 
peak in ultracold collisions of atoms and excited trimers 
in the nth state, since a^ g is quite well separated from 
other special values of a. Of course, in real experiments 
the resonance position may deviate from the universal 
value ^ due to finite-range effects; however, even with- 
out the dimer-trimer intersection, i.e., without the af^ 
point in Fig. [TJ the tetramer may still exist as it hap- 
pens in our model with n = 1. On the other hand, 
the observation of the resonant behavior in elastic or 
inelastic dimer-dimer collisions requires a very fine en- 
ergy resolution since the tetramer lies very close to the 
atom-trimer threshold. The situation is even more com- 
plicated if deeply bound dimer states are present such 
that the Efimov trimers become UBS with finite width 
and lifetime. According to Ref. (3( , the trimer width is of 
the order 77* 6„ with 77* being the atom-dimer inelasticity 



parameter; for example, 77* ~ 0.03 in the experiment of 
Ref. ^] performed at a > and 77* w 0.1 in the exper- 
iment of Ref. [13 performed at a < 0. Thus, the width 
of the trimer greatly exceeds the distance (i?^ Q — bn) be- 
tween the tetramer and the atom-trimer threshold, such 
that from the experimental point of view the nth trimer 
and the {n, 0)th tetramer are on top of each other. Fur- 
thermore, in the presence of deep dimers the atom-dimer 
scattering length Ad is complex with finite real and imagi- 
nary parts. For example, |Re^d/a| < 40 at 77* « 0.03 [3]. 
This should be sufhcient for the existence of the {n, 0)th 
tetramer at a slightly above a^ but excludes all the oth- 
ers (71,771 > 1). For larger 77* values even the (ti, 0)th 
tetramer may be absent. Thus, the experimental obser- 
vation of the considered tetramers would be extremely 
difficult. 

IV. SUMMARY 

We studied bosonic tetramers that have approximate 
dimer-atom-atom structure. They are unstable bound 
states in the continuum affecting the scattering pro- 
cesses in the four-boson system. Exact four-particle 
equations for the transition operators were solved in 
the momentum-space framework to describe the four- 
boson collisions. We accurately achieved the univer- 
sal limit by considering reactions involving highly ex- 
cited trimers but, for a given trimer, the results are 
restricted to the deepest tetramer only. We calculated 
the atom-trimer and dimer-dimer scattering observables 
to extract the tetramer position and width. In partic- 
ular, we determined the tetramer intersection with the 
atom-trimer threshold and demonstrated that in ultra- 
cold atom-trimer collisions it leads to a resonant enhance- 
ment of elastic and inelastic reactions. We studied also 
the reliability of the three-body dimer-atom-atom model 
for the atom-trimer collisions: it may be reasonable for 
the elastic scattering near threshold and the real part of 
the tetramer energy at \Ad\ >> a but fails for the inelas- 
tic reactions and the tetramer width. 
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